Import the ESRI shapefile of German districts and the district attribute table. Join the two data frames, transform the CRS to “EPSG:3035” and check your changes.
We want a first descriptive visual of the distribution of Covid-19 cases in Cologne and the surrounding districts. Calculate the number of Covid-19 cases in the last 7 days (cases_7days) by population (population) and multiply with 100,000 (in Germany usually called “7 Tages Inzidenzzahl). Call the column”covid7d_rate".
Select Cologne (district_id==5315), find the surrounding districts and plot only the Covid Rate in Cologne and the surrounding districts.
You can use the dplyr function bind_rows to cobine the information of cologne and the surrounding districts.
Save your data set of Cologne and its surrounding districts as a ESRI Shapefile.